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S  =  entropy 
t  -  time 
T  =  temperature 
u  =  velocity 

z  =  longitudinal  coordinate 


SUBSCRIPTS 


h  -  thermal 

1  ~  longitudinal 

o  -  zeroth  order  value 

p  :  prop-national 
■-  constant  pressure 

t  =  transverse 

T  =  constant  temperature 

v  =  viscous 

=  constant  volume 

-L  =  perpendicular 

||  =  parallel 

SUPERSCRIPTS 

o  =  special  case  for  k  =  0,  or  k  =  0 

r  z 

o  o 

^s-  =  non-dimensional  quantity 
vector  quantity 
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ABJ' 1';?  ACT 

The  Photo-acoustic  effect  occurs  when  a  chopped  b.  . i;n  of  li^ht  shin¬ 
ing  on  an  absorbing  gas  produces  sound  waves.  From  the  amplitude  of 
these  waves  one  can  infer  the  absorption  spectrum  of  a  known  sample,  or 
detect  trace  gases  in  an  unknown  sample.  While  the  signal  strength  is 
inversely  proportional  to  the  sample  chamber  volume,  in  some  cases  it  is 

advantageous  to  design  a  large  cell  as  an  acoustic  resonance  chamber,  and 

& 

modulate  the  excitation  beam  at  a  natural  frequency  of  the  chamber.  The 
resonant  amplification  factor  Q(ttf)  depends  upon  the  system  losses.  The 
dominating  loss  mechanisms  are  the  thermal  and  viscous  boundary  layer 
effects.  By  rescaling  the  local  coordinates,  one  can  include  these  los¬ 
ses,  neglecting  terms  smaller  than  the  boundary  layer  thickness,  simply 
by  modifying  the  boundary  conditions  for  the  non-dissipative  wave  equa- 
tion.  From  there  one  can  derive  an  equation  for  Q(hr)  applicable  at  all 
sample  chamber  resonances.  This  equation  will  be  useful  in  evaluating 
the  relative  advantages  of  resonant  and  non-resonant  photo-acoustic 
spectroscopy. 

An  experiment  conducted  with  a  Krypton/ Argon  laser  shining  on  a  one 
atmosphere  sample  of  dry  air  with  .$%  NO ^  enclosed  in  a  cylinder  resulted 
in  amplification  factors  less  than  those  predicted.  An  analysis  of  the 
data,  and  the  results  of  later  experiments  at  lower  pressures  indicated 
that  the  unpredicted  losses  were  probably  due  to  imperfect  cell  construc¬ 
tion. 
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Chapter  Is  INTdOUJGTiON 


fn  1830  Alexander  Graham  dell  dlucov.-red  the  photo-acoustic  effect 
whereby  light  energy  is  converted  into  co-ind.  He  found  that  intensity 
nodulated  light  focused  on  a  solid  substance  enclosed  in  a  sample  cell 
produced  audible  sound.  Later,  both  Tyndall  and  Hell  observed  similar 
results  shining  light  through  sample  cells  containing  gasses  and  liquids. 
The  light  energy  passing  through  a  rotating  light  chopper  is  periodically 
;orbed  by  the  irradiated  molecules,  and  rereleased  as  heat  causing  a 
periodic  pressure  pulse,  or  sound  wave.  The  energy  absorbed  is  greatest 
when  the  light  energy  (hy)  equals  one  of  the  energy  transitions  within 
the  irradiated  molecules. 

Fifty  years  after  Bell's  discovery,  Viengrov  used  the  photo-acoustic 

effect  to  construct  the  first  gas  concentration  measurement  instrument, 

and  to  record  the  first  absorption  spectra.  When  he  tuned  the  incident 

light  to  a  single  wavelength,  the  intensity  of  the  resulting  acoustic 

signal  was  proportional  to  the  concentration  of  molecules  in  the  sample 

with  that  corresponding  energy  transition.  When  he  scanned  the  light 

source  across  a  range  of  wavelengths,  the  intensity  of  the  acoustic  sig- 

13 

nals  corresponded  to  the  absorption  spectrum  of  the  sample. 

In  recent  years,  the  development  of  lasers  as  sources  of  intense 
wavelength  tunable  light  has  revived  interest  in  the  photo-acoustic  ef¬ 
fect.  Kerr  and  Atwood^  (1968)  conducted  the  first  experiments  with  a 
laser  illuminated  optoacoustic  detector  measuring  the  weak  atmospheric 
absorption  of  laser  radiation  by  water  vapor.  Kreuzer^  performed  the 
first  gas  concentration  measurement  with  a  laser  illuminated  optoacoustic 
detector.  The  increased  sensitivity  of  optoacoustic  detectors  due  to  the 
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use  of  lasers  has  made  them  useful  in  measuring  atmospheric  pollution 


concentrations  in  the  parts  per  billion  range,  and  gas  absorption  coef- 

-8  -1 

ficients  as  low  as  10  cm  .  One  may  also  use  the  optoacoustic  detector 
to  determine  the  energy  transfer  rates  between  vibrational  and  transla¬ 
tional  degrees  of  freedom  in  gas  molecules  by  monitoring  the  phase  lag 
between  the  modulated  light  and  the  pressure  variation. 

In  most  photo-acoustic  experiments  the  sample  cell  is  a  cylinder; 
the  laser  shines  along  the  axis  of  the  cell;  and  a  microphone  converts 
the  acoustic  signal  into  an  electric  signal.  The  resulting  signal  ampli¬ 
tude  is  proportional  to  the  gas  absorption  coefficient  and  the  incident 
light  intensity,  but  inversely  proportional  to  the  chamber  cross  section¬ 
al  area.  Sensitive  acoustic  cells  then  are  usually  as  small  in  diameter 
as  possible.  Practical  considerations  of  aligning  the  laser  beam  and 

accommodating  commercial  microphones  limits  typical  cell  diameters  to  25- 
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50mm.  However,  in  1973  Dewey,  Kamm,  and  Hackett  demonstrated  that  in 
some  cases  one  could  enhance  the  acoustic  signal  using  a  larger  diameter 
cell  designed  as  an  acoustic  resonance  chamber,  and  modulating  the  exci¬ 
tation  beam  at  one  of  the  natural  frequencies  of  the  chamber.  In  many 
cases  this  resonant  operation  is  better  than  non-resonant  operation.  Fbr 
example,  Shumate  et.  al.^  used  a  resonant  optoacoustic  detector  to  meas¬ 
ure  water  vapor  absorption  of  carbon  dioxide  laser  radiation.  They 
needed  a  large  volume  to  surface  area  ratio  cell  to  minimize  the  effects 
of  water  vapor  adsorption  at  the  walls. 

This  paper  provides  a  rigorous  treatment  of  the  problem  of  pressure 
signals  produced  by  a  distributed  heat  source  in  a  cylinder.  The  results 
should  be  useful  in  improving  the  interpretation  and  design  of  photo- 
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acoustic  experiments.  Specifically,  they  can  help  in  judging  the  rela¬ 
tive  merits  of  resonant  and  non-resonant  cells.  The  results  of  the  anal¬ 
ysis  are  general  equations  for  and  phase  shift  governing  all  the 

resonances  of  a  cylinder. 

In  most  current  photo-acoustic  spectroscopy  literature  the  emphasis 
is  on  laser  light  absorption  and  relaxation  mechanisms,  with  only  a  sim¬ 
ple  analysis  of  the  acoustics  portion  of  the  problem.  Often  the  acoustic 
equations  and  boundary  conditions  used  do  not  include  any  dissipation,  or 

the  system  losses  are  determined  experimentally  without  regard  to  their 
2.4. 5.9 

origin.  *  '  *  While  a  non-dissipative  analysis  is  adequate  for  predic¬ 
ting  the  approximate  resonance  locations,  the  results  give  physically 
unrealistic  infinite  responses  at  resonance.  Some  major  sources  of  sys¬ 
tem  dissipation  are; 

1.  Losses  due  to  the  compliance  of  chamber  walls 

2.  Dissipation  at  the  microphone  diaphragm 

3.  Losses  from  wave  scattering  at  obstructions 

4.  Free  space  viscous  and  thermal  losses 

5.  Viscous  and  thermal  boundary  layer  losses 

12 

6.  Relaxational  losses. 

A  sufficiently  rigid  cell  wall  material  can  make  reflection  losses 
due  to  the  compliant  chamber  walls  insignificant.  Dissipation  at  the 
microphone  diaphragm  is  similar  in  nature  to  that  due  to  flexible  walls. 
Consideration  of  microphone  losses  as  energy  dissipation  alone  however, 
is  insufficient  since  the  energy  lost  at  the  diaphragm  is  closely  linked 
with  the  ability  to  detect  the  pressure  waves.  For  a  microphone  with  a 

12 

small  surface  area  compared  with  that  of  the  cell,  the  loss  is  minimal. 
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The  losses  due  to  scattering  from  microphone  mounts  and  inlet  port:;  >-tc. , 
while  significant,  are  extremely  difficult  to  predict,  and  the  cell  de¬ 
sign  should  be  as  smooth  as  possible  to  minimize  them.  The  next  section 
will  show  that,  for  typical  coll  dimensions  and  pressures,  viscous  and 
thermal  losses  in  free  space  are  negligible  compared  to  visco is  and  ther¬ 
mal  losses  in  the  boundary  layer.  Relaxation  losses  occur  because  of  the 
finite  energy  transfer  rates  within  the  absorbing  molecules.  The  beam 
modulation  rate  must  be  much  slower  than  the  vibrational  relaxation  time 
of  the  gas  to  minimize  these  losses.  Accurate  predictions  of  dissipation 
due  to  relaxation  effects  depend  upon  precise  knowledge  of  the  energy 
transfer  processes,  and  are  not  included  in  this  paper. 

Chapter  2  demonstrates  that  one  can  include  the  effects  of  viscosi¬ 
ty,  thermal  conductivity,  and  compliant  walls  on  an  acoustic  wave  in  a 
problem  simply  by  modifying  the  boundary  conditions.  Chapter  3  derives 
these  boundary  conditions  for  a  closed  cylinder.  Chapter  4  solves  the 
resulting  eigenvalue  problem,  deriving  equations  for  the  resonant  frequen¬ 
cy,  the  Q,  and  phase  shift  of  each  cylinder  resonance.  Chapter  3  describes 
an  experiment  conducted  to  verify  those  equations. 
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Sviptor  Z:  •'•i'lA's'iONl'i  i.oviON^ 

The  alternate  compressions  aijil  ran  •factions  in  a  so  md  wave  corre¬ 
spond  to  hi  h  and  low  lemper.tture  areas.  In  a  ,;.u;  with  finite  thermal 
oon-luo*  ivity,  the  faster  hi  ;h  t-  i-rat  ire  mol*-..*  diff  ire  i:.fo  the 
cooler  r«  .--ions,  e.j’ial  iz5  n&  the  t«  .pern  tire  and  l*.-:  awning  the  amplitide  of 
the  .;ojnd  wave.  The  fluid  velocity  which  accompanies  an  acoustic  wave  is 
proportional  to  the  ;;r  aiient  of  pressure  or  temperature.  Viscosity  tends 
to  equalize  the  fluid  velocity  alone;  the  wave,  and  to  deintensify  the 
sound.  When  a  pressure  wave  strikes  a  wall  which  is  not  perfectly  rigid, 
the  surface  reflects  only  part  of  the  wave  energy,  and  absorbs  the  rest. 
The  ratio  of  the  pressure  to  the  .normal  fluid  velocity  at  a  point  on  a 
surface  is  the  acoustic  impedance  of  the  surface. 

One  can  describe  the  acoustic  behavior  of  a  fluid  as  the  sum  of 
three  waves;  a  propagational  wave,  a  thermal  heat  diffusion  wave,  and  a 
transverse  shear  wave.  The  thermal  and  transverse  waves  will  be  negligi¬ 
ble  outside  of  thin  surface  boundary  layers,  but  they  will  account  for 
the  majority  of  the  wave  power  lost  in  a  small  container  like  those  used 
in  photo-acoustic  spectroscopy.  The  thermal  and  shear  waves  appear  be¬ 
cause  the  propagational  wave,  although  it  dominates  the  acoustic  behavior, 
cannot,  by  itself,  satisfy  the  boundary  conditions  for  temperature,  tan¬ 
gential  velocity,  and  normal  velocity.  Since  the  viscous  and  thermal 
boundary  layers  are  very  thin,  one  can  consider  the  power  loss  to  occur 
at  the  boundaries,  and  can  include  the  compliant  boundary,  thermal  mode, 
and  shear  mode  effects  simply  by  modifying  the  boundary  conditions  on  the 
propagational  mode. 

For  the  small  disturbances  associated  with  acoustic  waves,  each  of 
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the  variables  describing  the  fluid  state  consists  of  an  equilibrium  part 
plus  a  small  time  varying  acoustic  part. 

Pressure:  P  +  p(x,t)  Temperature:  T  +  z(x,t) 

Density:  ^0+  S(x,t)  Entropy/unit  mass:  S  +  r  (x,  t) 

Velocity:  u(x,t) 


The  linear  equations  governing  the  wave  motion  are: 


L 

(O  at  1  P  V*  iA- 


0 


Conservation  of  mass 


,V 


:a 


(2)  p}t  ~  ~Vp ■+  uX)  -yoC\7 *  CVxit) 

—r 


Conservation 
of  momentum 


(3) 

(<0 

(5) 

where, 

r 

\ : 
* : 
Cf: 

V 

c\  = 


a  t 


K  vt 


^  -  djo  K6  ( 

-  -r-  k  <_  |0  / 


Continuity  for  heat  flow 

Equation  of  state 
Second  law  of  thermodynamics 


coefficient  of  viscosity 
bulk  modulus 
thermal  conductivity 
specific  heat  at  constant  pressure 
ratio  of  specific  heats,  Cp/Cy 


Pr\. 


Any  vector  function  of  space,  such  as  u\  has  two  parts:  a  longitudi¬ 
nal  part  u |  for  which  V  *  -0  »  and  a  rotational  or  transverse  part 

A  X 
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ut  for  which  V'Ut  O  •  Sir.ee  V  *  Cv7p}'0,  one  can  separate  the 
Navier-Stokes  momentum  equation  into  two  equations# 


(6)  /?  ^  :  \/.p  e >[  +  vLv'Ci^) 

(?)  ^^  =  ^Va  (  Vxafc) 


The  conservation  of  mass  equation  becomes 
(8) 

Equation  7  defines  u^,  which  is  the  shear  wave  mentioned  earlier.  The 
shear  wave  is  not  coupled  to  any  changes  in  fluid  pressure  or  temperature. 


One  can  use  the  conservation  of  mass  equation,  its  first  derivative 


with  respect  to  time, 

v-u  -  -J-U 

v  U/L  -  f  it 


X7 .  l-U’f  -  _  L 

v  btx  it*- 


and  the  definitions  of  §  and  cr  (equations  4  and  5)  to  reduce  the  system 
of  equations  above  to  two  equations  in  two  unknowns,  p  and 


Equation  11  expresses  the  longitudinal  velocity  in  terms  of  p  and  ~C> 


/°  It*  3  V7  (  f  '(4/0  q ( j° 


One  can  also  write  these  equations  in  terms  of  (characteristic  lengths, 


each  on  the  order  of  the  molecular  mean  free  path. 


(12)  X 


k  pC, 


=  A 

ir 


z  %f*A 


V  f 


A  cvt  = 


c  U 


V  )(. 


>-•  (-y 

^  £  v  C  (x  o  •jD 


0^UC  = 

r  u 


These  equations  represent  two  kinds  of  waves  depending  on  which 
equation  (13  or  14)  dominates.  One  can  find  the  equations  governing 
each  wave  by  considering  separable  simple  harmonic  motion, 

V  ~  "  ^  ft  ~  ~  llo  jO 


V  >p 

V"  r 


kS 


t  to  C 


and  solving  for  the  two  solutions  to  the  resulting  quadratic  in  k  , 

(16)  ( AijAA ^  0 


One  solution  corresponds  to  nearly  adiabatic  motion;  the  propagation  ve¬ 
locity  has  a  real  part  equal  to  c,  and  a  small  imaginary  part  correspon- 


8 


ding  to  the  energy  loss  terms  in  1  and  1^.  This  is  the  p ropagational 
mode,  and  it  is  governed  (to  first  order  in  1^  and  1^)  by  the  following 
equations: 


|  2 

V  p  -  -  k 

f 

ip 

k  r 

"  c  ck 

(17a) 

i  *  r 

k  ~  cZ  L 

i  •  ;f':  ( 

&  J 

(17b) 

if  ( 

/-rX)f 

0 

(17c) 

-  _  (  ,  L. 

—  \  c  t 

»  -  Sf.  )  V 
/ 

p. 

2 

The  other  value  of  k  is  large  and  complex,  indicating  rapid  attenuation. 

This  wave  corresponds  chiefly  to  heat  diffusion.  It  is  the  thermal  mode, 

and  is  governed  by: 

1  Z  ^  1 

V  Z  -  "k  Z  d-t 

(18a)  |< 

,  'l~-  (  uD 

'  ' 

(18b) 

f- 

(18c) 

oZY  f 
~  pc,  ^ 

:Kvr . 

The  plane  wave  solution  for  the  thermal  wave  is 

Ac  kx  -  i  u>t 

£ 
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bC  c 

The  thermal  boundary  layer  thickness  then  Is  on  the  order  of  / 

'/  oO 


(19) 


For  the  shear  layer 


A 


W 


u. 


=  A  (L 


i  kx  -  ( c )  t 


64. 
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a 

-  A  [  V  -zU  ^  -  0  X  -  t'u-it  ]  . 


The  viscous  boundary  layer  thickness  is  on  the  order  of  y- 


Kf 


(20) 


6^ 


2J^ 

CO 


Therefore,  while  the  propagational  mode  decays  as  exp(-/2.(  cj)  [" -^y  f  (^*'0.^70, 
the  thermal  and  shear  modes  decay  as  exp(-  ijjfih.  $  and  exP(“ 

ii  “6  -  o 

One  can  see  that  in  a  photo-acoustic  chamber,  —  =  0(10  ,  10  '),  that 

:  a 

the  thermal  and  shear  waves  exist  only  in  a  thin  area  near  the  boundaries, 
but  account  for  most  of  the  system  losses.  For  typical  photo-acoustic 
frequencies  and  chamber  sizes  the  attenuation  of  the  propagational  mode 
is  negligible,  and  one  can  assume  all  the  power  is  lost  in  the  boundary 
layer,  or  since  the  layer  is  so  thin,  at  the  boundary  itself.  Further 
analysis  will  neglect  terms  of  order  1^  and  retain  only  those  terms  of 
the  same  order  as  the  boundary  layer  thickness:  dy,  d^,  -jly* 
larger. 
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Chapter  3:  dO'JNDA <Y  CONJ1T iONJ 


Consider  the  cylindrical  conf ig iration  shown  in  Pigure  1 


x/  Pig.  1 

The  propagations!  mode  equations  are  now 


(21a) 


-  R(<~3  0(0)  eT“ 


(2  lb) 


cilO 

jo2- 


(21c) 


(21d) 


y:i  „ 


(21e) 


^  p  L  U)C 


If  the  boundary  surface  is  rigid  the  tangential  and  normal  velocities 
equal  zero  at  the  surface.  In  the  case  of  a  slightly  flexible  wall  one 
can  assume  the  surface  is  one  of  local  reaction,  with  acoustic  impedance 
^  to  account  for  a  nonnegligible  velocity  normal  to  the  wall.  Tangential 
velocity  would  still  be  negligible.  The  heat  capacity  of  the  boundary 
material  is  usually  much  greater  than  that  of  the  fluid  so  the  tempera¬ 
ture  variation  T approximately  equals  zero  at  the  surface.  Therefore,  a 


reasonable  approximation  to  the  true  boundary  conditions  is: 

(22)  U±  ~  U(l=  O  =  O 

at  r  =  a,  z  =  0,  z  =  L. 

But  the  propagational  wave  cannot  meet  all  three  of  these  conditions 
and  satisfy  its  own  equation  as  well.  One  must  add  enough  of  the  shear 
and  thermal  modes  so  that  the  combination  wave  satisfies  all  three  condi¬ 
tions.  A  discussion  below  of  the  r  =  a  boundary  will  demonstrate  how  the 
other  modes  affect  the  propagational  boundary  conditions. 

Assume  that  the  thermal  mode  has  the  separable  form 


=  Men  ^C9) 


-  LCjt 


(23) 


-  c  uJ> 


v  ^  =  4c 


-Zi 

z 
Vx 


The  substitution  r  =  a  -  d^x  rescales  the  local  coordinate  normal  to  the 

one  ti 

A'1  H 


boundary,  and  enables  one  to  eliminate  terms  smaller  than  d,  . 

h 


.IV 


■"  "  H 


(24) 


r- 


b 


Ci-.)' 


b 


'(' )n  As  r  ■ 


4  o  e 


)  <P  F 


Physically  the  solution  must  decay  with  positive  x;  so  b^  =  0.  The  ther¬ 
mal  solution  should  just  cancel  the  temperature  fluctuation  due  to  the 
pro pa  national  wave  at  r  =  a. 


Th  i  s  in.pl  :  <•;.  that 


-  t 


\> 
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•  >  rs~7  CI'l') 

(25)  A  ~  VZ 


■•'V, 


The  thermal  wave  has  with  it  a  corresponding  velocity  up  oc  Xp/r.  The 

2 

Q  and  z  components  are  of  order  d^,  but  the  r  velocity  component  is  not 
negligible. 

A  .  \  /"h  ~7  C.  I  -  *■  )  ~  jr~  "  £  ° ^ 


(26)  w^r  '  />C  6th 


7ik  Ruo  ( 1  -<■) &Z e 


The  shear  wave  will  cancel  the  tangential  velocity  due  to  the  propa- 
gational  wave. 

uV  rr 


(2?) 


VxVx  ut  =  Z  '  d, 


Letting  r  =  a  -  dyx  and  dropping  small  terms,  and  adjusting  constants  to 


cancel  the  $  and  z  components, 


_r.  J  .!  p  c'vry  p 
(28a)  =  top  73  /-  C 

(28b) 


y  (|- JC-J  S 


-  \  r  -  dv  / 

'  1  r  ■>  ^Cl'J 

U,  -■  t Kt K*  V  <r.J  ?  C 

t,  ^  , 


The  rotational  velocity  is  not  coupled  to  any  changes  in  temperature 
or  pressure,  but  it  does  have  a  component  in  the  r  directi  n  which  one 
finds  from  the  definition  of  u^, 

V  •  il  -  o 


(29)  U. 


^L)^'K)9z.d 


/  n! 


.  r  -o 
U-t)  d* 


—  1 >  { 


The  the iiii al ,  shear,  and  propagation;!!  components  of  the  normal  fluid 
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velocity  together  must  satisfy  the  u^  =  -p/^  boundary  condition 


(  u.,  •+  u 


r  + 


V 


)  =  -  fy? 


Solving  for  urp  results  in  the  true  boundary  conditions  for  p^. 


(30) 


U  r 


-X  _ 

^  a,  - 


i  ( 1  -  0 


>r 


1 


+  /°c 


+  kX'fcl: 


2c 


at  r  =  a 

The  three  terms  within  the  parentheses  correspond  to  compliant  boundary, 
thermal  conductivity,  and  viscosity  effects  respectively.  Assuming  the 
elastic  boundary  effects  are  small  (-^r  -  0[  ^w"] ),  equation  30  is  simply 


the  non-dissipative 


=  o 


a 


boundary  condition  with  small  modifications 


to  include  losses.  Since  thermal  losses  are  caused  by  compression  which 
is  related  to  the  normal  velocity,  the  thermal  effect  is  equivalent  to  an 
additional  acoustic  conductance.  The  viscous  term  is  more  complicated. 
Since  viscous  losses  are  related  to  the  tangential  velocity  of  the  fluid 
just  outside  the  boundary  layer,  the  viscous  term  involves  eigenvalues 
from  the  other  two  orthogonal  directions.  In  non-dimensional  form,  using 
the  conventions  below,  the  r  =  a  boundary  condition  is 


One  can  find  the  z  =  0  and  z  =  L  boundary  conditions  in  a  similar  manner 
using  the  substitutions  z  =  L  -  d^  yy  at  the  z  =  L  boundary,  and  z=  d^  ^y 
at  the  z  =  0  boundary.  The  z  =  0  and  z  =  L  boundary  conditions  are: 


(32)  r 


6\  L  L\) 


0-0 


at  z  =  0 


(33)  h- 


[/  (  cO  V  v 


\]/  y(  ,  Kr 


.  -v  ±j 

at  z  =  — . 
a 


The  boundary  conditions  couple  the  different  modes  of  vibration  to¬ 
gether.  The  eigenvalues  for  the  radial  modes  appear  in  the  longitudinal 
boundary  conditions,  while  the  longitudinal  and  azimuthal  eigenvalues 
appear  in  the  radial  boundary  conditions.  One  can  resolve  this  coupling 
problem  by  using  the  zeroth  order  approximation  to  the  radial  eigenval- 

V 

ues,  which  comes  from  the  =  0  condition,  in  the  longitudinal  boundary 
conditions  and  vice  versa.  This  works  because,  while  the  modified  bound¬ 
ary  conditions  produce  a  change  of  order  in  the  zeroth  order  eigen¬ 

values,  this  in  turn  produces  a  change  of  only  order  1^/a  in  the  other 
boundary  conditions. 


Chapter  4s  EIGENVALUES 


The  modified  boundary  conditions  (equations  31 • 32 , 33 )  can  be  applied 
to  the  nonhomogeneous  wave  equation  in  cylindrical,  coordinates  which  gov¬ 
erns  the  propagational  mode. 

V1-  =  FLc,i) 

While  the  forcing  function  K(r,t)  is  arbitrary  for  the  purpose  of  finding 

Q(io),  it  is  commonly  expressed  in  photo-acoustic  spectroscopy  literature 
1) 

as  -  ~~£  -tt’  ,  where  ^  is  the  net  rate  of  energy  release  into  the  ther- 
C  OL 

6 

mal  mode.  _ 

,  \  y  \  -  is 

W  [V1--  jp'f  -  <>  H 

One  begins  by  solving  the  homogeneous  problem  for  the  system  eigenfunc¬ 


tions. 

'V 

(35)  k  f  --  0 

The  asimithal  solution  is 


/%.  r 


k  =  kr  ■*  k 


t . 


©M  r  -*■  b„  >*\9. 


The  condition  ©,„  Co)  -  ©  (in)  imp)  ies  that  m  is  an  integer.  The  radial 

solution  R  is  a  linear  combination  of  Bessel's  functions  J  (k  r)  and 

m'  r 

^  v, 

N  (k  r).  The  coefficient  of  the  N  (k  r)  term  must  equal  zero  since  the 
m  r  m  r 

v 

problem  must  have  finite  solution:;  in  the  bounded  region,  and  N  (k  r)  is 


singilar  at  r  -  0, 

R 


C 


.....  ^  ^  k.  r  ) 


The  eig'Tiv.,1  i<  .  w".1  I  be  just  a  small  perturbation  of  order  iH^/a  from  the 


rigid,  nr>rivi:-.~o  nonco.vl  jetive  eigenv  1  ues . 


where  for  m  =  0 


k  =  0.000, 
r 

o 

for  m  =  1 

kr  =  1.842, 
o 

for  m  =  2 

\  =  3.054, 
o 

for  m  =  3 

k  =  4.201, 
ro 

etc. 

These  values  of  k 

and  z  =  — , 
a 


3.832,  7.016,  .  .  . 


5.332,  8.536,  .  .  . 


6.706,  9.970,  .  .  . 

8.015,  .  .  . 

r  appear  in  the  longitudinal  boundary  condition  at  z  = 
o 


0 


The 
are  small 


longitudinal  eigenvalue  k^  and  the  corresponding  eigenfunction 
perturbations  from  the  =  0  boundary  condition  values. 


Z 


k  -  where  p 

zo  L 

is  an  integer 

Applying  the 

boundary  condition, 

the  values  of  and 

X  to  first  order 

in  -<[lh/a  are 

(38a) 

''U  V 

2  B,  c  Z\  .) 

X  X 

(38  b) 

-  d>,  (Z-K.) 

X  ' 

whore, 

tV- 

■-  Ce;  (7  1(1 

\\  /  (  \  (  -4  .j^io  ] 

1  HO  Ov,.  1 

¥1). 

Both  //  an'l  5^ are  undefined  for  k  =  0.  The  solutions  for  k  =  0  are 

Zn  Z- 
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(39a) 


I 

! 


(39b) 

For  p  f  0 
(40a) 


t 


S' 


/  3  B2- 

—  i/^‘ 

CV  1/  « ^ 

S,(3„)  k  V/z 

2  J1 


Z-  rvi  n  i 


a 

U 


£Tva 


2-/)  i  rev.  1/ 


*2 


!-j  '■  l'  '  ' > ,,  /  k  2S1  k ) 

\  L_  L _  -p  ir  —1  k  _p  kb- 


./Ik 

i,'  6X 


For  p  =  0 


(40b) 


7 


»vl  *\0 


C&z. 


n  Vk  213,o,)ap„ 


f  T  ^ 

1 1  — 


L. 


The  k  =  values  appear  in  the  radial  boundary  conditions  at  r  =  a. 
Z_  Li 


At  r  =  1  , 


4? r  "iPu  b2(cj.)R 

csf  v~  I;  2.  - 


where 


’  K(Z)  -  CCJ  7  +  V  v  + 


8y  using  Taylor  series  expansions  around  k^  to  approximate  J^k^r)  and 

-  /7m \  at  k  +  k  -|  -1'  ,  one  can  show  that 
d(k  r )  r  1  >'  k- 

'  ■r  '  A 


fc 

V  2 

A/ 

■+" 

k 

r 

o 

to 

(41  a) 


L  [k ,  1 3 ) 

(tk-W) 


-v  2  2  ~ 

This  is  undefined  for  k  =  m  which  occurs  only  at  k  =  m  =  0. 

r  r 

^  o  o 

For  k  =  m  =  0 

r 

o 


(41  b) 


'  )  I  j  -A.  \ 

z.  J  Vt,;  )  !  /  > 


'  <r 


V  -  2  .  j 

'L  i 
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For  k  /  0 
r 

o 


(42a) 


-) 

\ 


k, 


^  i\  p 


0.  k,.‘  (V  z) 


r  j. 


? or  k  =  0 
r 

o 


(42b) 


R.„,p 


The  R,  and  2  eigenfunctions  correspond  to  a/.imithal,  radial,  and  lon¬ 
gitudinal  standing  waves  as  pictired  in  Figure  2. 


AZIMUTHAL 


m  =  0,1,2,.. 


RADI AL 


n  —  0,1,2,.. 

Fig.  2 


LONdITUDiNAL 

A 

h'  ^ 


f 


y 


A„„,,  e  1 


—  -  lO 


•nr-,  p  ""T  ) 

The  next  step  in  a  complete  solution  is  to  express  Q  in  terms  of  the 


eigenfunctions. 


*v»  n  t 


'  rV)  n  1  (  ^ 


However,  that  requires  knowledge  of  the  absorption  characteristics  of  the 
gas,  and  is  not  within  the  scope  of  this  paper.  One  can  get  some  very 
important  information  about  the  shape  of  the  resonant  peaks,  and  the 
phase  shifts  through  resonance  without  knowing  . 
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4^  * 


uf_.  _  -_LY-i)iCo_ 

w)  q  ■  5“-rk,N^') 

_  2  2  2 

System  resonances  occur  at  to  =  Re(k  +  k  ),  while  the  finite  imagi- 

r  z 

nary  values  of  f  and  keep  the  resonant  peaks  from  being  infinite. 

One  measure  of  the  sharpness  of  a  resonance  peak  is  the  Q(w)  at  reso¬ 
nance.  The  acoustic  Q(cj)  is  the  resonant  frequency  divided  by  the  peak 
width  at  the  half-power  point.  For  a  system  excited  by  a  single  laser 
pulse  the  rate  of  signal  decay  depends  on  the  magnitude  of  the  system 
dissipation.  Q(io)  =  2^  where  $  is  the  logarithmic  decrement  or  frac¬ 
tional  amplitude  loss  per  radian  of  the  signal.  The  acoustic  signal  amp¬ 
lification  due  to  resonance  reported  by  Dewey  et.  al.  is  proportional  to 

One  can  look  at  A  /b  near  resonance  by  defining  a  as  the  zeroth 
mnp'  mnp  °  * 

order  resonant  frequency  plus  a  small  perturbation  of  order  i/T^/a. 


:>  -  a,  -  -0  It 


s: V  t  v 


For  k  4  0,  k  /  0 
r  z 

o  o 

(45)  -A'-ijj  ...  -  -  . 

ia.,.p  (*/./ ok. k).k^.2  ^ 

In  most,  photo-acoustic  cells  the  walls  are  sufficiently  rigid  to  assume 

(_  A  0(l ,/a).  In  such  a  cane 

ft  t-T  ■  \ 


A . ,  . 

Tv . '  >0  U.’A 


| (.y~ o  t ) ,  i ((  .  i)<) -f)_i  tv 

v  I.  ,  ’  ,v  \j  1 V  .  .u  zj 


where 


Viscosity  and  thermal  conductivity  shift  the  cylinder  resonances  to 
higher  frequencies.  The  resonant  peak  occurs  at 


(4  7) 


LJ 


re  ■ 


=  4 


^  D 


To  find  Q (tj)  one  needs  to  solve  for -fly  where  I  /B  ,vi  ^  |I>  1  — 

'/n  .  To  first  order  in  1^/a  the  resonant  peak  is  a 

parabola  and 


(48) 


n.%  --  D  i  u>y  D. 


Therefore  the  width  at  the  half-power  point  is  2  »  and  the  formula 

for  Q(tj)  is 


(49) 


40re, 


cJr 


The  phase  angle  as  one  passes  through  resonance  is  equal  to  the  arc¬ 
tangent  of  the  imaginary  part  of  A  / id  divided  by  the  real  part.  Ihe 
-*  mnp7  mnp 

solution  for  the  phase  angle  0  is 

2  /x;  v.  y  ■:  s 

-//  x  -  X  -  Ut-  c^uU 


(50) 


^  -  £ 


AH 


For  the  special  cases  when  either  k  =  0  or  k  =0  only  the  formula  for 
t  rz 

o  o 

D  changes. 


r  or  k  /  0,  k  ~  0 

r  z 

o  o 


(51) 


n  -  i?  lx 


2> 

ro 


For  k  =  0,  k  /  0 
r  z  / 

o  o 


Some  literature  sources  refer  to  an  acoustic  Q(^j)  which  is  measured  at 
the  half-amplitude  rather  than  half-power  point  of  the  signal.  The  for¬ 
mula  for  the  half-amplitude  ^(a)  differs  from  the  half-power  Q(lj)  by  a 
factor  of  l/f3» 

Another  formula  for  Q(tu)  developed  from  Morse  and  Ingard^  for  the 

3,12 

m  =  0,  k  =0  radial  modes  is  common  in  the  photo -acoustic  literature. 
zo 

(53)  _ 1= _ _ _ 

Q^>  =  (/-O-V 14  k) 

This  formula  is  usually  adequate,  since  the  first  radial  mode  is  the  most 
strongly  excited  by  the  usual  beam,  cylinder  configuration.  When  equa¬ 
tion  49  is  put  into  the  same  form  as  53  for  easy  comparison  the  results 


"  |  p L.  k  Ck  <'<i)  V(j  a  I  (i-i  -i  -  ' 


O  W  k,  t  0 

ro 

1  /o.  k  =0 


k  r~ 

L 


p  TT  (X 


kr  -  the  r^gts  of  the  first  derivative  of 
&  the  m  order  Bessel's  function 


0 


Equation  'jk  reduces  to  53  for  m  =  k  = 

zo 

Equation  5^  is  more  general  because  it  is  valid  for  all  the  system 

resonances.  Current  research  requires  this  additional  flexibility.  For 

example,  some  multipass  configurations  which  bounce  the  laser  beam  back 

and  forth  inside  the  cell  to  increase  the  effective  absorption  cross 

section  may  not  excite  the  first  radial  mode  optimally.  Angus,  Marinero, 

and  Colles  observed  the  spectrum  of  NO ^  operating  in  a  small  diameter 

1 8 

cell  at  the  first  longitudinal  resonance.  They  achieved  an  amplifica¬ 
tion  factor  of  50  without  suffering  from  losses  due  to  a  larger  cell. 

Koch  and  Lahmann  detected  small  levels  of  sulphur  dioxide  in  air  using  a 
first  longitudinal  resonance.^ 


Chapter  5s  EXPERIMENT 


To  test  the  applicability  of  the  results  presented  in  the  first 
half  of  this  paper,  we  devised  an  experiment  to  measure  the  values  of 
Q(u))  and  phase  shift  through  the  acoustic  resonances  of  a  cylinder. 
Figure  3  summarizes  the  experimental  set  up. 


Fig.  3 

A  visible  laser  light  beam  was  passed  through  a  rotating  chopper 
wheel,  and  then  through  the  sample  cell.  The  resulting  pressure  pulses 
were  detected  by  a  microphone.  The  signal  from  the  microphone  was  am¬ 
plified  and  transmitted  to  a  lock-in  amplifier.  The  signal  was  synchro¬ 
nously  detected  by  the  lock-in  amplifier  which  was  referenced  to  the 
chopping  frequency.  During  the  experiment  the  chopping  frequency  was 
steadily  increased,  and  the  magnitude  and  phase  outputs  of  the  lock-in 
amplifier  were  plotted  on  the  chart  recorder  as  a  function  of  frequency. 
The  linearity  of  this  system  was  previously  tested,  and  it  is  assumed 
that  the  signal  plotted  was  proportional  to  the  pressure  changes  in  the 
cel  1 . 


F4 


We  conducted  our  experiments  primarily  at  atmospheric  pressure,  and 
room  temperature.  We  filled  our  cell  with  dry  air  and  nitric  oxide. 
Before  each  experiment  the  cell  was  pumped  out  overnight  with  a  vacuum 
pump.  We  filled  the  cell  with  1  to  3  Torr  of  NO,  monitoring  its  pres¬ 
sure  on  a  Bourdon  Gauge.  Then  we  filled  the  cell  to  one  atmosphere  with 
dry  air  by  passing  room  air  slowly  through  a  liquid  nitrogen  trap  to 
condense  out  all  water  vapor.  Upon  contact  with  air,  NO  reacts  to  form 

NO^  which  is  brown  in  color,  and  a  strong  absorber  of  visible  light. 

1 7 

Figure  4  shows  the  absorption  spectrum  of  NO^  •  and  the  locations  of  the 
output  line;  from  the  Spectra  Physics  model  171  Argon/Krypton  laser. 


Krypton  Lasers 


WAViiLENGTH  (A) 
Fig.  4 


We  used  both  the  Argon  and  Krypton  laser  tubes  during  different  parts  of 
the  experiment.  The  difference  in  wavelengths  between  the  two  lasers  made 
no  difference  in  the  experiment,  since  we  were  concerned  only  with  cre¬ 
ating  enough  signal  to  measure  Q(jj),  and  not  in  the  absorption  mecha¬ 
nism:;. 

2 

ted  a  Uruel  and  Kjaer  microphone  (sensitivity  =  4.93/dyne/cm  ) 


'We  mo  i  ri 


midway  along  the  axis,  flush  with  the  side  wall  of  the  cell.  A  ramp  vol¬ 
tage  applied  to  the  chopper  wheel  controller  increased  the  chopping  fre¬ 
quency  linearly  during  each  run. 

We  used  an  existing  acoustic  cell  (a  =  5.20?cm.  L  =  9.2075cm.)  with 
large  glass  windows  entirely  covering  the  ends  of  the  cell.  See  Figure  5* 


Such  a  conf ig  iratiori  made  it  difficult  to  distinguish  between  window 
absorption  and  gas  absorption  signals,  but  again  we  were  concerned  only 
with  having  enough  signal,  and  not  in  where  the  signal  came  from.  The 
large  windows  enabled  us  to  move  the  laser  beam  around  within  the  cell  to 
vario  isly  excite  and  suppress  different  modes.  Specif ically,  an  off  cen¬ 
ter  bean  could  excite  the  azimuthal  modes  which  were  not  excited  by  a 
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centered  beam.  Shining  the  beam  through  the  cell  at  r  =  .63a  suppressed 

the  first  radial  mode,  since  r  =  .63a  corresponds  to  the  first  zero  of 

the  J  Bessel's  function, 
o 

In  calculating  the  expected  Q(u)s  for  the  various  resonances  we 
assume  that  the  concentration  of  NO^  (<.5/0  was  so  small  that  the  dry  air 
used  as  a  buffer  gas  determined  the  cell  acoustics  and  dissipative  param¬ 
eters.  Table  1  lists  the  necessary  physical  constants  for  air  at  one 
atmosphere  from  the  Chemical  Rubber  Company  Handbook  of  Chemistry  and 
Physics. 

Table  1 


7V= 

1.40 

C  =  .241 

P 

Tempera ture  (C^) 

18.3 

21.1 

23.9 

26.? 

c  1 

(m/s) 

342.2 

344.0 

345.6 

347.3 

P  < 

kSm/lt)  „  z 

1.212 

1.200 

1.189 

1.178 

K  1 

'cal/gm  cm  (c °/cm)l0  ) 

60.81 

6l  .2? 

61.74 

62.20 

(micropoises ) 

182.0 

184.2 

185.5 

186.6 

Jfr 

.8514 

.8512 

.8509 

.8503 

JTh7a  1 

.081x1  O'5 

1 .088xl0"3 

1 .094x1 0"3 

1 .101x10 

Appendix  A  lists  the  mode  shapes  and  frequencies  of  the  30  lowest  natural 
frequencies. 

We  did  not  observe  all  of  the  resonances  listed  in  Appendix  A.  Many 
were  only  weakly  excited,  or  buried  tinder  the  strong  signal  of  a  nearby 
resonance.  The  odd  harmonics,  of  the  longitudinal  modes  were  especially 
difficult  to  detect.  Ihey  were  excited  by  the  attenuation  of  laser  power 
as  it  passed  through  the  absorbing  gas.  Kven  with  a  strong  absorber  like 
NO-,  the  odd  longitudinal  harmonics  were  much  weaker  than  the  radial  sig- 
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nals.  The  location  of  the  microphone  at  the  midway  point  along  the  axis, 


a  node  of  the  odd  longitudinal  harmonics,  made  detecting  them  even  more 
difficult. 

The  calculated  and  measured  values  of  Q(cj>)  for  the  observed  resonan¬ 
ces  at  one  atmosphere  are  in  Table  2. 


Table  2 


Mode 

Qcalc 

^Sneas 

Q  /Q 

hne  as  1 

(1.0,0) 

533.1 

383.4 

.719 

(2,0,0) 

5^9.2 

363.3 

.661 

(0,1,0) 

1135 

572.1 

.504 

(0,2,0) 

1543 

438.6 

.284 

(0,0,1) 

503.8 

412.1 

.818 

(0,0,2) 

714.2 

577.4 

.808 

(1.1.0) 

1279 

567.6 

.474 

The  signal  to  noise  ratio  for  this  experiment  is  probably  on  the  order  of 
100:1,  and  should  not  cause  any  error  in  measuring  Q(a>).  However,  there 
may  be  some  error  in  the  measuring  of  Q(gj)  due  to  uncertainty  in  fixing 
the  base  line  of  a  peak  from  which  we  determine  its  half  maximum.  This 
is  because  so  many  of  the  peaks  overlapped  with  peaks  of  other  resonances. 
Graphs  of  the  phase  and  amplitude  versus  frequency  for  each  resonance  are 


in  Appendix  B 


CONCLUSIONS 


The  results  for  the  longitudinal  modes  were  quite  close  to  their 
predicted  values,  while  the  azimuthal  and  radial  modes  were  not.  Some 
possible  sources  of  the  discrepancy  are: 

1.  Rough  Surfaces  -  With  boundary  layer  thicknesses  on  the  order  of 
_  3 

10  cm.,  slight  surface  roughness  could  have  seriously  affected  the 
results. 

2.  Compliant  Walls  -  While  the  cell  was  bored  from  a  heavy  block 
of  aluminum,  the  glass  windows  which  acted  as  end  walls  were  thinner  and 
less  rigid. 

3.  The  Microphone  -  The  microphone  diaphragm  was  23- 8  mm  in  diame¬ 
ter,  and  much  less  rigid  than  the  cell  walls.  The  mounting  of  the  flat 
microphone  was  not  flush  with  the  curved  walls,  which  could  have  signifi¬ 
cantly  altered  the  cell  geometry. 

4.  Volume  Losses  -  The  theory  in  chapter  2  predicts  these  losses 

_3 

will  be  on  the  order  of  10  times  the  boundary  layer  losses. 

Rough  surfaces  would  increase  viscous  losses.  The  windows  were 
smooth,  so  we  expected  to  see  rough  wall  losses  mostly  on  modes  tangen¬ 
tial  to  the  side  walls,  the  longitudinal  and  azimuthal  modes.  Since  the 
radial  modes  suffered  the  greatest  losses  beyond  those  predicted,  rough 
surface  losses  were  probably  not  the  major  source  of  error  in  this  exper¬ 
iment. 

Compliant  windows  would  show  up  as  additional  acoustic  conductance 
for  waves  propagating  normal  to  the  end  walls.  These  were  the  longitudi¬ 
nal  modes  whose  Q(l')s  were  very  close  to  their  predicted  values.  Compli¬ 
ant  windows  were  not  likely  to  be  the  major  source  of  error  in  this  ex- 


periment. 

Wave  reflection  off  the  flat  spot  on  the  side  wall  and  the  scat¬ 
tering  off  of  the  microphone  mounting  would  undoubtedly  have  the  most 
effect  on  the  radial  and  azimuthal  modes.  The  radial  and  azimuthal  modes 
have  the  largest  unexplained  losses,  so  the  microphone  was  probably  a 
major  source  of  experimental  error. 

One  way  to  verify  that  the  unpredicted  losses  are  due  to  imprecise 
cell  construction,  and  not  to  an  error  in  predicting  the  fluid  mechanical 
losses  is  to  conduct  the  experiment  at  lower  pressure.  This  raises  the 
relative  importance  of  the  fluid  mechanical  losses,  without  changing  the 
scattering  losses,  and  should  result  in  better  agreement  between  predic¬ 
ted  and  measured  Q(u>)s. 

We  measured  Q(oj)  for  the  first  radial  and  first  azimuthal  modes  at 

one-half  atmosphere,  and  Q  /Q  ,  improved  from  .504  to  .562  ,  and  .719 

iHG^S  G  a  1  G 

to  ,766  respectively.  The  small  improvement  was  very  close  to  that 
expected,  but  still  within  the  bounds  of  the  experimental  error.  We 
needed  further  experiments  at  l/4  and  l/8  atmosphere  over  a  wider  range 
of  resonances  to  confirm  our  hypothesis.  However,  the  laser  source  broke 
down  before  we  could  take  anymore  data. 

In  future  experiments  one  should  take  great  care  to  make  the  real 
cell  match  the  theoretical  geometry.  Mounting  the  microphone  flush  with 
an  end  wall  may  solve  the  scattering  problem.  The  experiment  should  cov¬ 
er  a  wide  range  of  pressures,  and  possibly  use  more  than  one  microphone, 
thus  providing  access  to  modes  with  a  node  at  one  microphone. 

This  paper  provides  a  needed  method  for  predicting  the  acoustic 
quality  of  a  given  resonant  photo- acoustic  spectroscopy  cell.  It  presents 


a  rigorous  treatment  of  the  resonant  responses  of  a  cylinder  to  a  dis¬ 


tributed  heat  source  by  showing  that  one  could  include  viscous,  thermal, 
and  compliant  boundary  losses  to  first  order  simply  by  modifying  the 
boundary  conditions  on  the  nondissipative  wave  equation.  The  results  of 
the  analysis  are  formulas  for  Q(gj)  and  phase  shift  through  resonance. 

This  paper  also  presents  the  results  of  an  experiment  conducted  to  test 
these  formulas.  Because  of  their  sensitivity  to  small  changes  from  the 
cylindrical  geometry,  the  formulas  failed  to  accurately  predict  the  exact 
response  of  the  system,  but  are  still  useful  for  comparing  the  relative 
advantages  of  various  experimental  designs  and  cell  configurations. 


APPENDIX  A 


The  Mode  Shapes  and  Frequencies  of  the  30  Lowest  Natural  Frequencies 
(r  =  5.207cm,  L  =  9.2075cm) 

m  =  the  azimuthal  index 
n  =  the  radial  index 
p  =  the  longitudinal  index 


(k2  +  k2)^ 

'  r  z 

(m,n,p) 

f  (18.3-26.7  C°) 

1.777 

(0.0.1) 

1860  -  1886 

1.842 

(1.0,0) 

1928  -  1955 

2.559 

(1.0,1) 

2678  -  2716 

3.054 

(2.0,0) 

3196  -  3242 

3-533 

(2,0,1) 

3698  -  3750 

3.553 

(0,0,2) 

3718  -  3772 

3.832 

(0,1,0) 

4010  -  4068 

4.002 

(1,0,2) 

4188  -  4248 

4.201 

(3,0,0) 

4397  -  4460 

4.224 

(0,1,1) 

4421  -  4484 

4.561 

(3.0,1) 

4773  -  4842 

4.685 

(2,0,2) 

4903  -  4973 

5.226 

(0,1,2) 

5469  -  5543 

5.317 

(4,0,0) 

556 5  -  5644 

5.330 

(0,0,3) 

5578  -  5658 

5.332 

(1,1,0) 

55 80  -  5660 

5.502 

(3,0,2) 

5758  -  5841 

5.606 

(4,0,1) 

5867  -  5951 

5.639 

(1,0,3) 

5902  -  5986 

6.143 

(2,0,3) 

6429  -  6521 

6.395 

(4,0,2) 

6693  -  6789 

6.40? 

(1,1,2) 

6705  -  6801 

6.416 

(5,0,0) 

6715  -  6811 

6.564 

(o,l,3) 

b870  -  6968 

6.658 

(5,0,1) 

6968  -  7068 

6.706 

(2,1,0) 

7018  -  7119 

6.786 

(3,0,3) 

7102  -  72 04 

6.937 

(2,1,1) 

7260  -  7364 

7.016 

(0,2,0) 

7343  -  7448 

7.106 

(0,0,4) 

7437  -  7543 

J. 


APPENDIX  B 


GRAPHS  OF  PHASE  AND  AMPLITUDE 
VERSUS  FREQUENCY  FOR  MEASURED  RESONANCES 

The  amplitudes  of  all  the  peaks  are  normalized  to  one. 
The  frequency  is  given  in  Hz. 
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